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In addition to the inclusive cross sections discussed within the QCD-parton model, in the 
regime of multiple parton interactions, different and more exclusive cross sections become 
experimentally viable and may be suitably measured. Indeed, in its study of double parton 
collisions, the quantity measured by CDF was an "exclusive" rather than an inclusive cross 
section. The non perturbative input to the "exclusive" cross sections is different with respect 
to the non perturbative input of the inclusive cross sections and involves correlation terms of 
fN| , the hadron structure already at the level of single parton collisions. The matter is discussed 

$_i ' in details keeping explicitly into account the effects of double and of triple parton collisions. 
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INTRODUCTION 



The growing importance of multiple parton interactions (MPI) at high energy has stimulated a 
lot of interest in the phenomenon, in view of the forthcoming results at the LHC(lM2|L MPI are es- 
sential to describe the features of the minimum bias and of the underlying event 3f] 0] jl] p | f^] [8] [9] 
and may represent an important background in many channels of interest at the LHC[1][2], not 
only in processes where the cross section is large, as the unitarity issue at the origin of the effect 
might suggest, but also in cases where the cross section is rather small, like for the search of the 
Higgs boson 10 1 U| or in the case of the production of equal sign W boson pairs 12|[13|[2|. On 



o , 

Q\ , the other hand MPI are by themselves an interesting topic of research, since by studying MPI 

one may obtain informations on the multi-parton structure of the hadron [HQ]. 

Up to now the direct observation of MPI has not been easy. The direct measurement of 
^ ■ MPI requires in fact the identification of the final fragments of the multiple process and the 

reduction of statistics, due to the request of large momentum exchange in each hard interaction, 
has restricted considerably the possibilities of a direct study of the phenomenon. To measure 
the MPI one needs moreover to separate the background due to hard radiation. A given multi- 
partons final state may in fact be produced either by a multiple or by a single parton collision. 
The separation between the two contributions has proven to be experimentally feasible in the 



case of double parton collisions [14] 15] [16] [17(] - The enhanced contribution of MPI at high energy 
will facilitate considerably direct studies of the phenomenon and one may reasonably expect that 
the separation of the two different contributions will be done more easily at the LHC, at least in 
the simplest cases of MPI. 
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In the regime where MPI may be observed directly, in addition to the inclusive cross sections 
usually considered in large momentum exchange processes, one has the possibility to measure di- 
verse and more exclusive cross sections, computable in perturbation theory and linked differently 
to the hadron structure jl8||. In pp interactions MPI are dominated by independent collisions, 



initiated by different pairs of partonsjjjj] 20]. A direct consequence is that the multi-parton 
inclusive cross sections are proportional to the moments of the distribution in the number of 
collisions On the other hand, a statistical distribution may be characterized either by its 
moments or by its different terms. While the moments of the distribution in the number of 
collisions are measured by the multi-parton inclusive cross sections, the different terms of the 
distribution are measured by a different set of observables, which one may call "exclusive" cross 
sections. Inclusive and "exclusive" cross sections result from independent measurements and are 
linked in a different way to the hadron structure. The two sets of cross sections are however con- 
nected by sum rules. By testing the sum rules, namely by looking at the number of terms needed 
to saturate the sum rules in a given phase space region, one measures the effects of unitarity 
corrections, which allows to control the consistency of the analysis and provides an additional 
handle to obtain information on the multi-parton correlations of the hadron structure. 

Interestingly, in its study of MPI, the CDF experiment did not measure the inclusive cross 
section of double parton scattering. The events selected where i fact only those which contained 
just double parton collisions, while all events with triple scatterings (about 17% of the sample of 
all events with double parton scatterings) where removed[l7j]. The resulting quantity measured 
by CDF is hence different with respect to the inclusive cross sections usually discussed in lar ge m 
physics. In fact it represents precisely one of the "exclusive" cross sections recently discussed [181] - 

While the inclusive cross sections are linked directly to the multi-parton structure of the 
hadron the link of the "exclusive" cross sections with the hadron structure is much more elabo- 
rate. The requirement of having only events with a given number of hard collisions implies that 
the corresponding cross section (being proportional to the probability of not having any further 
hard interaction) depends on the whole series of multiple hard collisions. The number of hard 
partonic collisions which can be observed directly is nevertheless limited, which allows to discuss 
the "exclusive" cross sections by expanding in the number of elementary interactions. Simple 
and directly testable connections between inclusive and "exclusive" cross sections may hence be 
established by saturating the sum rules, which link "inclusive" and exclusive cross sections, with 
a finite number of terms. 

The purpose of the present paper is to discuss the "exclusive" multiparton scattering cross 
sections, going up to the third order in the number of collisions and keeping two-body parton 
correlations explicitly into account, while the effects of the three-body parton correlations in the 
hadron structure will be neglected. Explicit expressions in terms of the two-body correlation 
parameters will be derived in a few simplest cases, where correlations will be assumed to depend 
only on the transverse coordinates. 



2. "EXCLUSIVE" CROSS SECTIONS 



In pp collisions, the inclusive cross sections are basically the moments of the distribution 
of the number of MPI[18(]. The most basic information on the distribution in the number of 
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collisions, the average number, is hence given by the single scattering inclusive cross section of 
the QCD parton model. Analogously the i^-parton scattering inclusive cross section cr^-gives the 
Kih moment of the distribution in the number of collisions and is related directly to the -ftT-body 
parton distribution of the hadron structure: 



a K = J D A (xi...x K ]bi...b K )cr(xix 1 ) ...o{x K x' K ) 

xDb(x'i . . . x' K ; b\ — (3 ... bx — (3)dxidx[d 2 bi . . . dxKdx' K d 2 bKd 2 f3 
{N(N - 1) . . . (N - K + 1)) 



K\ 



-Chard (1) 



where ahard represents the contribution to the total inelastic cross section due to all events with 
al least one hard interaction, while D(x\ . . . xk', b\ . . . bx) is the K-partons density of the hadron 
structure, with transverse parton coordinates b\ . . .bx and fractional momenta x\ . . . xk, (3 is 
the hadronic impact parameter and a the parton-parton cross section. 

A way alternative to the set of moments, to provide the whole information of the distribution, 
is represented by the set of the different terms of the probability distribution of multiple collisions. 
Correspondingly, in addition to the set of the inclusive cross sections ajc, one may consider the 
set of the "exclusive" cross sections <7/v, where one selects the events where only N collisions are 
present. One hence has: 



O^hard = a N ' ° K = IF] N ^ ' 

N=l N=K 

which represents also a set of sum rules connecting the inclusive and the "exclusive" cross sections. 

While the non perturbative input to the inclusive cross section o~k is given by the iT-parton 
distributions of the hadron structure, as implicit in Eq.(2), the non-perturbative input to the 
"exclusive" cross sections is given by an infinite set of multi-parton distributions. The request 
of being in a perturbative regime limits however the number of partonic collisions and the sum 
rules in Eq.(2) are saturated by a few terms, in such a way that the "exclusive" cross sections 
can be expressed by finite combinations of inclusive cross sections. A particular case is when all 
correlations C n with n > 2 are negligible (higher order correlation terms may be introduced in 
the picture of the interaction as sketched in the appendix). In that instance all sums in Eq.(2) 



can be performed[2l||[22| and (neglecting all rescatterings) the hard cross section (Thard can be 



expressed in the following functional form: 



0~hard(f3) 



1 



expj— J dudu'dja(u,u')dj>^ Za[J\Zb[J' 



(3) 



j=J'=i 



where u = {x, 6}, v! = {x' , b' — f3} and a(u, u') is the interaction probability of the two partons 
with coordinates u and u'. For simplicity flavor indices are omitted and the integration limits are 
set by the limits of the phase space window where the MPI are observed. The two-body parton 
correlations are introduced through the functional 
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Z[J + 1] = ex p|y D(u)J(u)du + ~J C{u,v)J{u)J{v)dudv\ 
= ^2^\J J i u i) ■ ■ ■ J ( u n)D n (ui . . .u n )dui . . .du n 



(4) 



which generates the non-perturbative input of the n-parton inclusive cross sections, the inclusive 
n-body parton distributions Dn(u-\ . . . u„). Under these conditions Eq.(3) can be worked out 
fully explicitly. One obtains 2 11] [221] 



VhardiP) = 1 - exp -- ^a n - ^y^&n/n 



1 



L 2 



(5) 



where 



«n = (-l)^ 1 ^ DA(ui)a(u 1 ,u[)CB(u[,U2)a(u2,u 2 )CA(u2,u 3 ) . . . 

n 

■ ■■&(u n ,u' n )DB(u r n )Y[duidu' i 



(6) 



i=l 



and the chain, which starts with A, may end either with A or with B, depending wether n is 
odd or even. For b n one has 



b n = J C^UnjU^&fa^^CBiu'x,!^) ■ ■ ■ 

n 



(7) 



1=1 



and in this case only even values of n are possible. 

The exponential in Eq.(5) represents the probability of no interaction at a given impact 
parameter j3. All "exclusive" cross sections can be obtained from the argument of the exponential. 
One may start from the partonic interaction probability 



i-IId 



(8) 



where in the product each index assumes a given value only once, in such a way that possible 
re-interactions are not included. The probability of having only a single interaction is expressed 

by 







T 9 \ na-^) 



g=l kl ij^kl 



(9) 



3=1 



while the probabilities of a double and of a triple collision are 



1 ( d\ 2 n 1 n 

2! \~~d~q) n ( j - ^ = ^ n ( j - 

V y / i,j=l 9 =i ' kl rs ij^kl,rs 



1 _ d_ 

3! I % 



n (i-^v) 



9=1 



9=1 



1 n 



3! 



fci rs tu 



ij^kl,rs,tu 



9=1 



and the corresponding expressions for the "exclusive" cross sections 



UK™ 



9=1 



9=1 



9=1 



^{[X'( 9 )f-X"(g)}e- X (9) 
^{X'"( 9 ) + [*'( 5 )] 3 - 3X'(g)X"(g)} e ~ X ^) 



9=1 



where X = \{^ a n + Yl K/n). 

It's convenient to expand X and its derivatives in the number of elementary collisions 

X = X 1 +X 2 + X 3 + ... 



where 



X, 



Xo 



J Djy(u)a(u,u')DB(u f )dudu' 

2 

D j4 (ui)(t(ui,u / 1 )C , b(u / i, u' 2 )a(u 2 , u 2 )Da{u 2 ) JJ duidu^ + A <-> i? 

i=i 

1 /■ 2 

~2 / C '^(' u i'' u 2)o-(«i,'"i)C'B(«i,'"2)o-('"2,«2) JJduidu^ 



i=i 



X 3 = I D A (u 1 )a(u 1 ,u 1 )CB(u' 1 ,u' 2 )a(u 2 ,u 2 )C A (u 2 ,u 3 )a(u 3 ,u 3 ) 
3 

xD B (u 3 ) Y[duidu'i 
i=i 



The derivatives at g = 1 give 
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X[(v,v') 
X' 2 {v,v') 

X' 3 (v,v') 



D A (v)a{v,v')D B (v') 

— Da(v)ct(v,v') J CB{v\u' 1 )d(u' l ,ui)DA{ui)duidu' l + A <-> B 

— j Ca(ui,v)ct(v, v')Cb(v' , Ui)a(u'i,ui)duidui 
D A (v)a(v,v') / Cb(v',u' 1 )&{u' 1 ,ui)Ca(ui,U2)&{u2,v! 2 ) 



xDb(u' 2 ) Y\ diiidu'i 
i=l 

+ J DA(ui)a(ui,u' l )CB(u' l ,v')a(v',v)CA(v,U2)v(u2,u 2 ) 

2 

xDb(u' 2 ) Y^duidu'i 
i=i 

+ J D A (u 1 )a(u 1 ,u 1 )C B (u' 1 ,u 2 )a(u 2 ,U2)CA(u2,v)a(v,v') 

2 

xD B (v') Y^diiidu'i 



i=i 



(14) 



and 



X 2 (v l ,v' 1 ;v 2 ,v 2 ) = - DA{vi)a(v 1 ,v' 1 )CB(v' 1 ,V2)a(v 2 ,V2)D A (v2) + A ^ B 

-C A (v2,v 1 )a(v 1 ,v' 1 )C B (v' 1 ,v 2 )a(v 2 ,V2) 
X' 3 \v 1 ,v > 1 ;v 2 ,v 2 ) = 2 DA(v 1 )a(v 1 ,v' 1 )C B (v' 1 ,v' 2 )a(v 2 ,V2) J C A (v 2 , ui)a(ui, u[) 

xDB(u'i)duidui 
+Da{vi)ct(vi,v' 1 ) J CB(v[,u' 1 )a(u' 1 ,ui)CA(ui,V2)duidu' 1 

xa(v2,v' 2 )DB{v2) 
+ J D A (ui)a(u 1 ,u f 1 )C B (u f 1 ,v' 1 )du 1 du l a(v' l ,vi)CA(vi,V2) 

xa(v 2 ,v 2 )D B (v 2 ) 

^3 , {v 1 ,v' 1 ;v2, v' 2 ; v 3 , v 3 ) = 6D A (vi)a(vi, ^)C B (^, i/ 2 )a(v 2 , v 2 )Ca{v2, v 3 )a(v 3 , v' 3 )D B {v 3 ) 

(15) 

By substituting the expansions in the number of elementary collisions in the expressions of the 
interaction probabilities and by expanding the exponential, one obtains the expressions: 
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a[ = (X[+X' 2 + X' 3 )(l-X 1 -X 2 + X 1 -X 1 /2) 
2xa 2 = (X[ ■ X[ + 2X[ ■ X' 2 - X 2 - Xg){l - Xi) 

3 x 3 " = \{X' 3 " + X[.X[.X[- 3X[ ■ X 2 ') 



(16) 



where a[ etc. are the "exclusive" cross sections, differentiated according with (14) and (15). The 
integrated "exclusive" cross sections hence are 



<7i = X\ — X\ — X\X 2 + Xf/2 + 2X 2 — 2X 2 X\ + 3X 3 

2 

3 x 03 = 3X 3 + (Xi) 3 /2 - 3X ± X 2 



2xa 2 = Xf + 4X X X 2 - 2X 2 - 6X 3 - Xf + 2X ± X 2 



(17) 

The sum rules of Eq.(2) are satisfied as follows 

01 + 2 X 2 + 3 x 3 = Xi 

-Xf + 2X 2 + Xf - 2X 2 - X X X 2 + Xf/2 - 2X 2 X X + 3X 3 
+4X X X 2 - 6X 3 - Xf + 2X1X2 + 3X 3 + (Xi) 3 /2 - 3XiX2 
= Xi = as 

2 x 02 + 6 x 03 = Xf- 2X 2 

+4XiX 2 - 6X3 - Xf + 2X1X2 + 6X3 + (Xi) 3 - 6X1X2 
= Xi 2 - 2X 2 = 2 x D 

6 x 03 = 6X3 + (Xi) 3 - 6X1X2 = 3! x T (18) 

where 05, a B and ar are respectively the single, double and triple parton scattering inclusive 
cross sections. Explicitly 



05 = Xi = J D A aD B 

a D = ^[Xf - 2X 2 ] = ^[J D A aD B ■ D A aD B + J D A aC B aD A 

+ / D B aC A aD B + / C A aC B a 



= X - J [D A D A + C A ]aa[D B D B + C B ] (19) 

where [DD + C] = D 2 , the two body parton distribution as defined in Eq.(4). An analogous 
expression may be written for ot- 
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The relations (18) may be inverted 



5\ = CTg — 2d £) -\- 3&T 

a 3 = ot (20) 

which allow to express the scale parameters characterizing the double and triple parton collisions 
in terms of the single scattering inclusive cross section as and of the single and double parton 
"exclusive" cross sections 5\ and o 2 ■ 



1 0% 



2(J eff 

1/ .1,1 

3 6 ra 2 eff 

where the scale factor of the triple parton scattering cross section has been characterized by the 
dimensionless parameter r. 



3. CORRELATIONS IN TRANSVERSE SPACE 

Multiple parton collisions are most important in the region of small fractional momenta, 
where the large population of partons may dilute correlations and e.g. correlations due to energy 
conservation may not be of major importance. A simplest possibility is hence to neglect alto- 
gether correlations in fractional momenta and to work out some case where only correlations in 
transverse space are present and which allows an analytic treatment. To disentangle the effect of 
correlations in transverse space from other sources of correlation, we will consider the instance 
where all correlation terms give zero, when integrated on the transverse coordinates, in such 
a way that all other variables remain uncorrelated and, in particular, the distribution in the 
number of partons is Poissonian. 

The actual dependence of the correlation on the transverse variables is not prescribed by 
general principles, in this section two choices are presented, a Gaussian shape and an exponential 



shape, both were used in a previous discussion[18|. 



3.1. Gaussian density 



Using Gaussian distributions for the parton densities in transverse space and for the correla- 
tions one obtains closed analytic expressions: 
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D{x,b) = G(x)f(b) 

f(b) = g(b,R 2 ) 

C(x 1 ,x 2 ;b 1 ,b 2 ) = G(x 1 )G(x 2 )h(b 1 ,b 2 ) 

Hh,b 2 ) = c-g(B,R 2 /2)h(b,X 2 ) 



h(b,x 2 ) = ±g{b,\ 2 h) 



7=1 



(22) 



where G{x) represents the usual one-body parton distribution, g(b, X /j) = rj ■ g(b, X /j) 



g(b,R 2 



exp(-b 2 /R 2 ^ 



TV. 



R 2 



and 



B = [bi + b 2 ]/2 

b = [bi - b 2 ] (23) 
in such a way that the following relations hold 

J d 2 bg(b,R 2 ) = l, J d 2 b 2 g(b 1 -b 2 ,R 2 1 )g(h 2 ,R 2 )=g(h 1 ,R 2 + R 2 ) 

J Hb^b^dH = (24) 

One may define the correlation length r c as the value of b where the correlation term h(b,X 2 ) 
changes sign. With our definition of h(b, A 2 ) one has r c = X. 

To define unambiguously the "correlation strength" c, one needs to normalize properly the cor- 
relation term h(b, X 2 ) (which integrates to zero). Our choice is 

f h(b,X 2 )d 2 b = l 

J\b\<r c 

which gives r] = e, where e the Euler's number. 

The integrations on the transverse variables of the terms in Eq.'s 13-15 give 
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D A oD B - J d 2 bd 2 (3g(h - (3, R 2 A )g(h, R 2 B ) = 1 
D A &D B ■ D A aD B -► y dVhcPpgib! - (3, R A )g{h u R 2 B )g(b 2 - (3, R 2 B )g(h 2 , R 2 B ) 



D A aD B ■ D A aD B ■ D A oD B ^ d 2 b x d 2 b 2 d 2 b 3 d 2 {3g(h x - (3, R 2 A )g(h u R 2 B )g(h 2 - (3, R 2 A ) 

xg(h 2 ,R 2 B )g(b 3 - (3, R 2 A )g(h 3 , R 2 B ) 



-^{R\ + R B Y 

D A aC B aD A - ^ d 2 b x d 2 b 2 d 2 {3g{h x - (3, R A )h B {h u h 2 )g{h 2 - (3, R 2 A ) 

_ c B e X B 

= —(2R 2 A + \ 2 B ) 2 

D A aC B aD A ■ D A oD B ^ d 2 b x d 2 b 2 d 2 b 3 d 2 /fy(bi - /3, i&)Mbi, b 2 Mb 2 - A #4) 

x 5 (b 3 -/3,i?i) 5 (b 3 ,i?|) 
_ c B e 2A| 



3vr 2 (^ + J R|)(2i?2 +A 2)2 
C A aC B <7 y d 2 b 1 d 2 b 2 d 2 f3h A {h 1 - (3,b 2 ~ P)h B (b u b 2 ) 

= c A c B e 2 2X 2 A X 2 B 
vr (A 2 + A|)3 

C A aC B a ■ D A aD B -► ^ ' d 2 b 1 d 2 b 2 d 2 b 3 d 2 : /?Mbi - /3,b 2 - /3)/i B (b l5 b 2 ) 

x 5 (b 3 -/3,i?i) 5 (b 3 ,i?|) 
_ c^cge 2 4A^A| 

3vr 2 (#2+^2+^)3 

D A aC B aC A a D B - ^ d%d 2 b 2 d 2 b 3 d 2 Pg^, i?i)Mbi -(3,h 2 ~ (3)h A (h 2 ,h 3 ) 

xg(b 3 - (3,R 2 B ) 
16c A c B e 2 2 2 2 
= 3^P|^' SB ' r) (25) 

where = (X/R)a,b, r = Ra/Rb and 



,2 2 

77/2 2 2\ _ ^ 



32r 2 + 100 + 32r~ 2 + 6s|(l + 4r~ 2 ) + 34(34 + 2(4r 2 + 1)) 



{4[34 + 2(4r 2 + 1)] + 2 [6 + 4(1 + 4r~ 2 )] 



Using Eq.s(20) and (21), all inclusive and "exclusive" cross sections, up to the triple order in the 
number of parton collisions, are expressed in terms of the single scattering inclusive cross section 
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as, of the "effective" cross section a e ff and of the parameter r. In the case of collisions of two 
identical hadrons the explicit expressions of a e ff and r are 



and 



16 x 3s 2 



(7, 



eff 



8itR 2 



1 + c • e , ~ ' ~" + c 2 • e 2 



(4 + 3s 2 ) 2 



35" 



(26) 



TO 



eff 



16vr 2 ii 4 



f 16 X 9S 2 no 

|l + c-e - , „_ ONO +c 2 -e 2 



(4 + 3s 2 ) 2 
rrelation 1( 

measured in the generalized parton distributions 23] 24]. 



2 ^/3s 2 3s 2 

—a + 6 x 64£4 , ,1 

3s 2 V 2 ' 2 ' 



(27) 



where s = X/R (A = r c is the correlation length) and i? 2 = |i? 2 is the square hadron radius 



3.2. Exponential density 

The case where the parton density has an exponential shape is better displayed in Fourier- 
transform representation: 



g(b,R 2 ) 



1 



(27T) 5 



-ikb 



1 



[l + k 2 R 2 ) 2 



d 2 k 



(28) 



while all other terms defined in Eq(22) are redefined according with this unique change. In 
particular, the correlation term is 



h(bi,h) = c- g(B,R 2 /2)h(b,X 2 ) 



ikb 



(kXf 



(l + k 2 X 2 ) 3 



d 2 k 



and g(b, X 2 /^/) = r\ ■ gib, X 2 /j). As in the previous case, one defines the correlation length r c as 
the value of the distance |t»i — h>2 1 where the correlation h(b, A 2 ) changes sign. The relation with 
the parameter A is r c = xoX with xq ~ 2.387 18y. The normalization parameter rj, defined by the 
requirement 



L 

J\b\<r c 



h(b,X 2 )d 2 b = l 



is now rj ~ 3.456. 

The integrations on the transverse variables cannot be always displayed in closed form, a 
relevant simplification is got by taking the parameters R and A to be equal in A and B. 
In this case the integrations on the transverse variables of the terms in Eq.'s 13-15 give 
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D A aD B - y d 2 bd 2 f3g{h - 0)g(b) = 1 



D A aD B ■ D A aD B 



d 2 b x d 2 b 2 d 2 (3g(h x - /3)<7(bi)<7(b 2 - 0)g{h 2 ) 
1 1 



4vri? 2 7 



D A aC B aD A - y d 2 b 1 d 2 b 2 d 2 pg(b 1 - /3)^(bi,b 2 ) 5 (b 2 - /3) 



4vri? 2 



7]F 



3 + 44a - 36a 2 - 12a 3 + a 4 + 12a(2 + 3a)ln(a) 
F(a) = 3(a-l)6a 

C A &C B o -► y d 2 b 1 d 2 b 2 d 2 f3h(h 1 - /3, b 2 - /3)/t(bi, b 2 ) 

c 2 2r/ 2 



D A aD B ■ D A aD B ■ D A aD B 



AirR 2 15s 2 

_L 1_ /"°° _(1 + x + y)[(l + x + y) 2 + 6xy]dx dy 

(4tt) 2 Wo (1 

1 



(4vri? 2 ) 2 



+ x) 4 (l + y) 4 [(l + x + y) 2 - 4xy] 7 / 2 
x 0.030 



D A aD B ■ D A aC B aD A 



'(4 



r 2 _ 

vr) 2 « 4 Jo (1 



(1 + x/4 + y) 2 - xy/2 



+ x/2 + y) 3 [(l + x + 2y + (x/2 - y) 2 ] 3 / 2 
dx dy 



s 2 y 



C A aC B a ■ D A aD B 



(4ttR 2 ) 2VH 



2 2 _l±i r 

77 4vr 2 i? 4 s 2 J 



(l + s 2 y) 3 (l + x) 4 (l + x/2) 2 
1 



dx 



y 2 dy 



(l+x) 4 (l+x/2) 4 (l+y)6 



c 2 r? 2 

~ 2 X - 0256 



(4irR 2 Y s 



1 /' 7r d(j) f°° dx 
D A aC B aC A aD B -► cV / 4— / — — — — 

(4vri? 2 ) 2 7 vr J (1 + x) 2 (l - 



dy 



w 



w 



(1 + u>) 3 (1 + w') 



rj L 



(4itR 2 ) 



+ x/2) 2 (l+y) 2 (l+y/2) 2 



(29) 



where s = X/R, w = s 2 [x + y/4 — ^/xycoscf)] and it/ = s 2 [x/4 + y — y / xy cos 0] . 

The functions F, if and L defined above, are given in Table I for four different values of the 
parameter s 2 . 

The root mean square radius R in now given by R 2 = 12R 2 , so a more meaningful reference 
parameter is s = r c /R = xq/\/Y2 ■ X/R. As a function of R, s and of the correlation strength c, 
the expressions of the effective cross section and of the parameter r are 
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Table I: The functions F, H and L in Eq.'s(29) for three different values of s 2 



2 


0.5 


1 


2 


3 


s 2 


0.2372 


0.4744 


0.9488 


1.4232 


F[l/s 2 } 


0.0140 


0.0667 


0.2735 


0.5834 


H[l/s 2 ] 


0.0104 


0.0120 


0.0119 


0.0110 


L[l/s 2 } 


0.0089 


0.0094 


0.0078 


0.0064 




TO 



eff 



0,27 
^R4 



l + c-mr,H(^^)+c 2 -r] 2 



xl 0, 64 



(31) 



3.3. Superposition of Poissonians 



A particular case, where all n-body correlations are important and which can be worked out 
explicitly, is when the parton distribution is given by the superposition of different Poissonians. 
The superposition of Poissonians is naturally obtained when introducing diffraction in a multi- 
channel eikonal model of high energy hadronic interactions 25] 26] 27]. To have some indication 
on this case we consider the simplest possibility where the probability P n , to find n partons 
within a given kinematical range, is given by the sum of two Poissonians with average numbers 
ni and ni 



Pn 



"1,2 



n 

7— r e 



' ~ n l + (1 - 7)^ e ~ n 2 
n! 



"1,2(2;, b)dxd 2 b, "1,2(2;, b) = G(x)g(b, R\ 2 ) 



(32) 



here 7 gives the relative weight of the two Poissonians and the integration limits in x are defined 
by the kinematical range relevant to the case of interest. Notice that, as we want to disentangle 
the effect of correlations in the transverse coordinates, while ni(x,b) 7^ n2(x,b), the integrated 
values n\ and rii are equal. The average density of partons with fractional momentum x and 
transverse coordinate b is 



00 „ 
(n) = ' S ^nP n = \jni + (1 — 7)n 2 ] = / D(x,b)dxd 2 

n=l J 

D[x,b) = GW[ 75 (^?) + (l-7)fl(^2)] 



(33) 



while for the average density of pairs of partons with coordinates x%, b\ and X2, 62 one obtains 
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oo 

(n(n - 1)) = n ( n ~ ^ Pn = b n ? + ( X - 7 ) n 2] 

n=l 

= J D 2 (xi,bi;x 2 ,b 2 )dxidx 2 d 2 bid 2 b 2 ; 
D 2 (x 1 ,b 1 ;x 2 ,b 2 ) = G(x 1 )G(x 2 )[7g(b 1 ,Rl)g(b 2 ,Rl) + (l- 1 )g(b 1 ,I^)g(b2,R 2 2 )] (34) 
and analogously 

£> n (xi,&i . . .x n ,6 n ) = 6r(»i) . . . G(x n ) 

x [70(61, fl?) • • • g(b n , R\) + (1 - 7)5(^1, Rl) ■ ■ ■ g(K, Rl)] (35) 

The expression of the inclusive cross section of N independent parton collisions apj hece is 

a N = ±a»x{ 7 2 J[g(P,2R 2 )} N d 2 (3 + 2 7 (1 - 7) f [g(fl , R\ + Rl)] N d 2 [5 

+(1- 7 ) 2 | [g((3,2R 2 )] N d 2 p} (36) 

The actual calculations are carried out for the case of the Gaussian parton density, where the 
effective cross section and the parameter r of the triple scattering inclusive cross section are 
given by 

iir^-^ 2l{1 -' ,) ^ +(l - 1? 'w^ (37) 

and 

4 -am • + ™ -7) ■ © 2 + a - ^ • p., 

where we made the positions 

i?f = a/? 2 , i?l = (2 - a)R 2 and R 2 = ^R 2 

with ^ 2 the mean square hadron radius measured in the generalized parton distributions. 

One recognizes that the three different terms in the curly brackets are the contributions to the 
scale factors due to all possible combinations of the different sizes of the two interacting hadrons 
(Ri-Ru R1-R2+R2-R1 and R 2 -R 2 )^, 

To make contact with the general formalism previously discussed, one may identify the cor- 
relation term by the relation 



D 2 (xi,b 1 ;x2,b 2 ) = [D(x%,bi)D(x 2 , b 2 ) + C(x\, 61; a? 2 , 62)] 



(39) 
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which gives 

C{xu h;x 2 , b 2 ) = G(x 1 )G(x 2 ){ [ ig {b u R 2 )g{b 2 ,R 2 ) + (1 - l)g(bi, R 2 2 )g(b 2 , R 2 2 )] 

-[ 7 5(6i,fi?) + (i-7M&i, J R 2 2 )] 

x[ 7 g(b 2 ,Rl) + (l-^g(b 2 ,Rl)]] 

= 7 (1 - j)G(x 1 )G(x 2 ){ [g(h,R 2 ) - g(h,R 2 )] 

x[g(b 2 ,R 2 )-g(b 2 ,R 2 )]] (40) 

The correlation strength c is hence expressed as a function of the relative weight of the two 
Poissoninas 7 by the relation 



c = 7 (l-7) (41) 

while the correlation length r c , now defined by the change of sign of the correlation terms for b\ 
or b 2 = r c , is given here below as a function of the mean square root hadron radius R and of the 
parameter a, which controls the relative value of the two transverse radii R\ and R 2 : 



R 



r< 2 



a)ln 



a 



(42) 



27l | corresponds to c 



\ and r r 



ln3 ~ .h2R. With R2 = .42fm' 



28] 



The case discussed in 
one obtains cr e // ~ 30mb, too large to explain the value of a e it observed by CDF 

Notice that in the two-channel eikonal model discussed in[25j] one obtains a value of a e jf in 
agreement with the experimental indication. The reason is that, in the model, compact hadronic 
configurations are characterized by a stronger Pomeron coupling, which corresponds associating 
a higher partonic population to the compact configurations. In the present case the distribution 
in the number of partons is, on the contrary, the same in the two configurations with transverse 
distances R\ and R 2 (namely n\ = n 2 , after integrating on b). One may hence conclude that the 
experiment indicates that the fluctuation of the whole hadron structure in its transverse size, 
to the extent suggested by diffraction, is not enough to explain the value of <7 e // ; which may, 
on the contrary, require the introduction of correlation terms of the kind discussed here above 
in subsections 3.1 and 3.2, or/and of correlations in fractional momenta. In this last instance 
the multiparton distribution of the hadron structure is different from a Poissonian, also after 
integrating on the parton's kinematical variables. 



4. CONCLUDING REMARKS 



Multiple parton interactions are going to play an important role at the LHC, both in the 
description of the properties of the minimum bias and of the underlying event and as a background 
to various channels of interest for the search of new physics. The study of MPI represents 
moreover the basic handle to obtain information on unknown non-perturbative features of the 
hadron structure, namely the correlations between partons. 
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In the present paper we have tried to identify the quantities which are most suitable to obtain 
the information on the non-perturbative features associated to the presence of MPI, in the case 
where the MPI are dominated by independent collisions initiated by different pairs of partons. To 
this purpose, in addition to the inclusive cross sections considered till now in hard processes, we 
made use also of the information provided by the "exclusive" cross sections. Following a previous 
article where the "exclusive" cross sections were introduced[l8(|, we have hence analyzed MPI 
considering systematically all terms up to triple scatterings. The inclusive and the exclusive 
MPI cross sections are linked by the sum rules of Eq.(2) and, by checking the number of terms 
needed to saturate the sum rules, one has a direct control on the importance of the unitarity 
corrections. The case where no more that three parton collisions give significant contributions 
leads to very simple relations between inclusive and "exclusive" cross sections, the relations (18) 
and (20) which, being a consequence of the definitions of the cross sections, hold rather in general. 
One may thus obtain the values of the relevant parameters, a e ff for the double parton collisions 
and the dimensionless parameter r for the triple (having defined in that case the scale factor as 

A convenient way to measure the scale factors may be through Eq.s(21), which make use only 
of single and double collisions terms. Notice that if the sum rules of Eq.(2) are saturated with 
three terms in a given phase space window (and hence in a given interval of x values) Eq.s(21) 
must hold. By measuring as, a% and 02 as a function of x in the given interval, one may hence 
obtain a reliable information on the dependence of the correlation terms on x. 

In our approach correlations are introduced in the most general way, as deviations of the multi- 
parton distributions from the Poissonian. Rather than trying to propose definite correlation 
models, our philosophy is hence to identify the observable quantities which are most suitable to 
obtain information on the correlation terms of the hadron structure. To have an idea of where 
the correlation parameters (correlation length and correlation strength) may be most relevant we 
have considered a few simplest cases. Of course correlations will depend on all variables and in 
particular on fractional momenta, because of conservation laws. Nevertheless conservation laws 
may not play a very important role when the parton population is large, namely at small x. In 
Section 3 he have worked out in full detail three simplified instances, where the dependence of 
correlations on x may be neglected and which allows a full (or almost full) analytic treatment 
(Gaussian and exponential parton densities and correlations, multi-parton distribtuions given by 
a superposition of Poissonian). The relevant non perturbative information, namely the quantities 
a e ff and r, are given as a function of the correlation parameters in Eq.s (26, 27), Eq.s. (30, 31) 
and Eq.s (37, 38) in the three cases. Notice that the hypothesis of a negligible dependence of 
correlations on fractional momenta is easily tested experimentally by looking at the dependence 
of a e ff and r on x. 

In our discussion we did not allow for the differences between partons (gluons and quarks, 
different flavors, valence and sea). When considering definite reaction channels, the relations 
obtained have hence to be adapted, taking into account that the information on correlations will 
have to be related to the different kinds of initial state partons involved in the interactions. 

APPENDIX 

When looking at the effects of the correlation terms of the original multi-parton distribution 
in processes where three or more partons undergo hard scattering, a natural question is how to 
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deal with higher order correlation terms. In this appendix we sketch the procedure to deal with 
this problem. When taking in consideration three-body correlations the starting point is: 



Z[J + 1] = exp j J D 



(u)J(u)du + — I ('(a. r). /((/)./( i')<lii(h' 



H — / T(u,v,w)J(u)J(v)J(w)dudvdw 
6 



y^^j y J(ui)...J(u n )D n (u 1 ...u n )du 1 ... 



du r . 



The possibility of giving a close form to Z does not exist any more, but a general procedure 
is available, and well known 29]: every functional admitting a formal series expansion can be 
expressed as 



<&[J] = $[6/6x] exp / J{u)x{u)du 



so the previous expression may be rewritten in the following form: 

Z\J + 1] = exp - f T(u, v, w)— -rf— ■ / . dudvdw 
L6 J o(u) d(v) o(w) 

xexpjy '[D(u) + x( u )]J( u )du + i y C(u, w) J(u) J(u )(indf;l 



The second exponential will give rise to terms similar to the ones in Eq. (6,7), the only difference 
is in the terms a n , which contain D + % instead of D. The whole expression is suitable for an 
expansion in T. When acting with the derivatives one finds that every term T is connected either 
with three terms M(u,v), defined by the relation below, ending in turn on a density D or with 
two of them, one of which ends on one D and the other is closed on the same term T. 

a n = (-l) n+1 J D A (u 1 )a(u 1 ,u' 1 )C B (u' 1 ,u' 2 )a(u2,U2)CA(u2,U3) . . . 

it 

• • • a(u n ,u' n )D B (u' n ) Y[ didn't 
i=l 

/n 
£>A(«i)^(tti,Mi)M n (ui,\)(T(n n , u' n )D B {u' n ) duidu^ 

i=l 

In particular the first two terms are: 

T(u, v, w)a(u, u')D(u')a(v, v')D(v')a{w, w')D(w')dudvdwdu' dv' dw' 
and J T(u,v,w)a(u,u)a(v,v')C(u',v')a(w,'w')D(w')dudvdwdu'dv'dw' ; 
and one expects the first to be the most important. 
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